Running coupling and screening in the (s)QGP 
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Emphasizing the importance of renormalization in the context of thermal field theory in general, it 
is pointed out that the Debye mass in the hot quark gluon plasma is determined by the coupling at 
the scale mo, not T as commonly presumed. The mended result agrees quantitatively with lattice 
QCD calculations in the strong coupling regime almost down to T c . 
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I. INTRODUCTION 

It is common knowledge that the presence of a medium 
does not introduce new UV divergencies in quantum field 
theoretical calculations, basically because the momenta 
of the medium fluctuations are constrained by the tem- 
perature T and/or the chemical potential /i. At times 
this fact is used to comfort a somewhat pragmatic strat- 
egy in practical calculations, namely discussing only the 
matter contribution to the quantity under consideration, 
deferring the renormalization (of the vacuum part) as 
standard 'textbook physics'. Such a procedure, how- 
ever, should be applied with the appropriate care in cases 
where new momentum scales occur, as in the context of 
statistical field theory. A related aspect thereof - which 
is crucial for phenomenology - is the arising ambiguity in 
the coupling, whose momentum scale then is not fixed by 
the calculation, but has to be chosen somewhat ad hoc. 

A widely used expression obtained from such a prag- 
matic prescription is 



<D 



= (l + in/)47ra(Q|)T 2 



(1) 



for the Debye mass in a hot quark gluon plasma (QGP) 
with rif light quark flavors. Describing the screening of 
(chromo-) electric fields, the Debye mass is a key param- 
eter for a variety of phenomcnologically relevant quan- 
tities. For quantitative estimates, the coupling a has 
prevalently been fixed at a 'typical' thermal scale of the 
order of the (bosonic) Matsubara frequency, Qt = c-2wT, 
with c ~ 1. Recently, the expression (1) has also been 
utilized as a parameterization of lattice QCD (1QCD) re- 
sults [1, 2]. It there turned out that the Debye mass 
in the non-perturbative regime near the (phase-) tran- 
sition temperature T c can be reasonably described by a 
re-scaling Ansatz nrhr), see below. 

A reckoned explanation of the observed enhancement, 
K ~ 1.5, is the next-to-leading order correction [3] to 
the result (1) which, under reasonable assumptions for 
the magnetic mass, is positive. As an aside it is noted 
that such an interpretation is quite non-trivial. Naively, 
one could expect that if mjj/T ~ {Aita) 1 ' 2 = g changes 
notably in a certain temperature range and the next- 
to-leading order correction is of a similar size as the 
leading term, one cannot parameterize the modification 
by a constant factor. However, in QCD the correction 
is logarithmically enhanced due to an IR sensitivity of 



the Debye mass, implying n nl ° = 1 + Ag (ln g 1 + c) , 
where A = 3/ (4tt(1 + \n f f/ 2 ) [3]. Thus, from A K nl ° = 
(dK nlo /dg)Ag = (n nl ° - 1 - Ag)Ag/g, a sizable change 
in g can indeed be reconciled with an almost constant 
enhancement factor if K nl ° — 1 — Ag w 0. Although this 
condition is fulfilled in the cases studied below, it turned 
out that the next-to-leading order correction can only 
partly account for the observed excess [2]. 

Pursued here is a more elementary interpretation of 
the large Debye mass near T c , namely by questioning the 
(temperature dependence of the) coupling in (1). In order 
to first define rigorously the notion of the running cou- 
pling for the following considerations, some basic facts 
about the coupling renormalization (at T = 0) are recol- 
lected in Sec. II. Then in Sec. Ill, the pole of the propa- 
gator at T > is expressed in terms of the renormalizcd 
coupling. By specifying the general result for the case of 
the Debye mass it is shown that the presupposed scale 
Qt in the coupling in Eq. (1) is not justified. In some 
places the notation is kept schematic, both for simplicity 
and to indicate that the considerations are, in principle, 
not restricted to QCD (albeit consequences for weakly 
coupled systems are only marginal). Although the frame- 
work is entirely perturbativc, it will be instructive - with 
regard to the so-called strongly coupled QGP (sQGP) 
- to confront the results with 1QCD calculations. The 
paper concludes with a summary and some implications. 



II. RUNNING COUPLING 

To relate the bare coupling a to an experimental 
observable, consider a scattering process in a kine- 
matic regime dominated by the t-channel contribution. 
Considering for simplicity a massless theory, the rele- 
vant part of the matrix element is a/ (P 2 — II vac (P 2 )). 
To 1-loop order and in dimensional regularization, the 
boson selfenergy has the generic form II vac (P 2 ) = 



ln(— P 2 //^ 2 )). Then, for a specific P 2 



aP 2 (e 
the matrix clement reads 



_ a(t e ) 



t e 1 - a (e- 1 - ]n(-i e /p 2 )) t e 



t, 



(2) 



The right hand side can be considered as experimental in- 
put, introducing the coupling a(t e ) at the scale t e , which 



2 



is related to the (infinite) bare coupling by 
a- 1 (t e ) = a' 1 - er 1 + ln(-i e //i 2 ) 



(3) 



An equivalent relation holds for the coupling ait) at an 
arbitrary scale i, consequently 



L (i e ) + 47r/3 ln(t/*e 



Here the previously omitted prefactor has been re- 
introduced, which determines the asymptotic behavior 
of the specific theory under consideration. For QCD 
with N c = 3 colors and nt light quark flavors, the lead- 
ing order coefficient of the (3- function is /3q = (^jr N c — 
|n/)/(47r) 2 > 0. The momentum dependence of the 
'running' coupling is fully specified by its value at the 
scale t e , or in a common alternative way of writing, 



*(*) = 



1 



1 



(4) 



4^/3o ln(|t|/A2) ' 

by the QCD parameter A. It is also noted that the matrix 
element (2) can now be expressed simply as 

a _ a(P 2 ) 

p 2 -n vac (P 2 )~ p 2 ' 

with an obvious diagrammatic interpretation as depicted 
in the insert of Fig. 1. 

Before turning to the case of finite temperature, it is 
instructive to challenge the 1-loop expression for the run- 
ning coupling by a comparison to non-perturbative 1QCD 
results. To that end, Eq. (4) is applied in the position 
space by replacing \t\ — > r~ 2 . On the lattice, the cou- 
pling can be obtained from the heavy quark potential 
V(r) or more conveniently, as advocated in [1] and refer- 
ences therein, from the corresponding force, 



,(r) 



3 r2 dV(r) 
4 dr 



The comparison in Fig. 1 shows an agreement which is 
at least as good as that obtained in [1] from the Cornell 
parameterization of the potential, V cp (r) = — | a/r + ar, 
with a fitted constant a and the string tension ^fa = 
420 MeV. Worth to note is that the adjusted value A w 
205 MeV for the running coupling (4) is verily in the 
expected ball park. This demonstrates that the pro- 
nounced increase of a qq (r) is not an unambiguous in- 
dication of non-perturbative effects, as a cp (r) — > for 2 , 
but can be utterly accommodated within the asymptotic 
form, at least at intermediate distances. This is a mo- 
tivation that phenomenological considerations based on 
the 1-loop running coupling may be viable even for mod- 
erate interaction strength, say a 0.6 (corresponding to 
r ^ 0.3 fm) as relevant in the following. 



III. DEBYE MASS 

The Debye screening mass can be calculated from the 
position of a pole of the dressed propagator in the pres- 
ence of a thermal medium. Here the selfcnergy receives 




FIG. 1: Running coupling at T = 0, for ly =2. Compared 
are the 1QCD results [1] (given here without error bars) to 
Eq. (4) with t — » r -2 (full line); the hatched band depicts a 
variation AA = 40 MeV. The dashed (black) line represents 
the result [1] from the Cornell parameterization of V(r), which 
for r J; 0.3 fm almost coincides with the full (red) line. 



an additional finite contribution, 

II = a [P 2 (e- 1 - ln(-P 2 /^ 2 )) + f(p ,p)] , (5) 



which in general depends separately on po and p. The 



pole is determined by p$ 
also be written as 



P 



IL(po,p) = 0, which can 



a 



ln(-P 2 /p 2 ) + f(p ,p)/P 2 



Recognizing in this expression the running coupling 
a(P 2 ), cf. Eq. (3), this relation is cast in the intuitive 
form 



P 2 = a(P 2 )f(p ,p). 



(6) 



While simple, this equation is remarkable since it is man- 
ifestly finite. The pole is determined only by the matter 
part of the selfenergy, whereas the (divergent) vacuum 
contribution is 'absorbed' completely in the renormal- 
ized coupling. The fact that the pole depends only on 
the measurable parameter A (while the auxiliary scale n 
has dropped out) is a strong formal indication that the 
approximation (6) has a direct physical meaning [11]. 

As a special case of Eq. (6), the Debye screening mass 
is defined as the pole position of the longitudinal bo- 
son propagator Don for Po = 0, at space-like momentum 
— P 2 = m 2 D . Then it is first obvious that screening is 
an effect solely due to the thermal fluctuations. The vac- 
uum contributions of the selfenergy (5) do not have to be 
'omitted', e.g., by misleadingly arguing that for large T 
they are exceeded (after renormalization) by the matter 
part, / ~ T 2 3> p 2 ■ On the contrary, they are actually 
necessary to render finite the resulting equation. Not sur- 
prising, in fact, is that the relevant scale in the running 
coupling is given by the Debye mass itself, viz 



m% = |A C (1 + \n f ) 4na(m 2 D )T 2 



(7) 
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Here the same perturbative result for f(0,p) as in Eq. (1) 
has been used, which for small p is momentum indepen- 
dent and explicitly gauge invariant [5]. Thus already to 
leading order (i. e. without using a resummed selfenergy 
calculated from dressed propagators), due to renormal- 
ization the Debye mass is determined by an implicit equa- 
tion. With the coupling (4), the self-consistent solution 
can be given in terms of Lambert's VF-function [12], 



bT 2 



W(bT 2 /A 2 ) ' 



(8) 



where b — ^(1 + gn/)//3o- Even without having to refer 
to the tabulated properties of the IT-function, two basic 
features of the solution can be easily inferred directly 
from Eq. (7). First, mu/T is a monotonically decreasing 
function of T/A and second, the self-consistent solution 
is larger than the expression (1) if mo < Qt- 

Although also Eq. (8) is justified a priori only for 
T ^> A, it is remunerating to compare it to 1QCD re- 
sults. In Monte Carlo calculations, the Debye mass can 
be obtained by two methods: from analyzing the large- 
distance behavior of (i) the heavy quark free energy F 
or (ii) the static electric propagator D 00 . For quenched 
QCD as well as for pure SU(2), the Debye mass has 
been calculated within both approaches. Even though 
for SU(3) the results are not fully consistent, see Fig. 2, 
they are clearly underestimated, by the aforementioned 
factor k ~ 1.5, by the expression (1) supplemented with 
the 2-loop running coupling [13]. The new formula (8), 
on the other hand, can describe the 1QCD results almost 
down to T c by appropriate choices of the dimensionless 
parameter A = A/T c . I will focus on the data sets [1] 
since they, in the quenched limit, match nicely with the 
findings in SU(2) [6], cf. Fig. 3. In this computationally 
less demanding case the lattice results derived from the 
free energy and from the propagator are indeed compat- 
ible. On the other hand, for n; = the Debye mass (8) 
is also iV c -independent. 

Adjusting the scale parameter in Eq. (8) to the QCD 
(JV C = 3) lattice results (excluding a narrow interval 
above T c ) yields A( n/=0 ) ~ 1 for the quenched case and 
\n /= 2) ~ 1-3 for two light flavors. One might suspect 
that the lattice results could be accommodated within 
this strictly perturbative approach only at the expense of 
the interpretability of the parameter. I thus emphasize 
that the fitted A are in a remarkable agreement with the 
known ratio of the QCD scale to the transition temper- 
ature [1], A/T c = 1/1.14(4) and A/T c = 1/0.77(21) for 
n/ = and rif = 2, respectively. To round up the emerg- 
ing consistent picture it is underlined that A w 205 MeV, 
as adjusted in Sec. II for the n/ = 2 QCD potential at 
T = 0, is (together with the commonly expected value of 
T c as 160MeV) also compatible with A(„ /=2 ). 

Due to the monotonic behavior of mrj/T alluded to 
before, Eq. (8) cannot reproduce the characteristic de- 
crease of the Debye mass very close to T c . A possible 
interpretation has been given as a Boltzmann suppres- 
sion of the time- like excitations (quasiparticles), which 
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FIG. 2: The QCD Debye mass for n } = and n, = 2. The 
1QCD results labeled F and Doo are from [1] and [2], respec- 
tively. The hatched bands depict Eq. (1) with the 2-loop cou- 
pling fixed at Qt = c • 27rT with c £ [|, 2] and values of A/T c 
as specified in the text. The dash-dotted (black) lines repre- 
sent the re-scaling fits [1] with the Ansatz mho; here c = 1. 
The thicker (red) lines show the self-consistent result (8) for 
some representative choices of the parameter A = A/T c . 
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FIG. 3: Comparison of quenched 1QCD results: SU(3) [1] vs. 
SU(2) [6]. The deviations very near T c are not unexpected 
due to the different order of the phase transition. 



become heavy near T c [7]. Within the present approach, 
it could alternatively indicate a non-monotonous a(Q). 

For large temperatures, the deviations of the results 
from Eqs. (1) and (8) decrease logarithmically, reaching 
the order of 10% only at 10 2 T C , cf. Fig. 4. From the 
slopes of the curves it can be expected that the re-scaling 
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FIG. 4: Top: Large-T behavior of the Debye mass for nf = 2 
(thin lines) and rif = 3 (assuming A( nf=3 ) = A( n , =2 ); bold 
lines). Bottom: The adjusted 1-loop coupling (4) at the scale 
mu vs. the 2-loop coupling at 27rT. 

Ansatz nrhu will overestimate the correct Debye mass at 
high temperatures (a hint thereof may in fact be spotted 
already in Fig. 2). Shown also in Fig. 4 is the prediction 
for rif = 3 massless flavors, which should constrain the 
Debye mass in the physical case (rif = 2 + 1). 

IV. CONCLUSIONS 

In a wider sense, it was demonstrated that there is no 
coupling in a (thermal) medium per se - it always is to 



be related to a specific quantity/scale. The correct scale, 
which should be an outcome of the calculation rather 
than 'chosen', is crucial in particular for larger coupling. 
The relevant scale for the Debye mass, 




is tojj itself, not the typical loop momentum Qt ~ 2ttT 
as often presumed. It's the interplay between vacuum 
and thermal fluctuations (accounted for by renormaliza- 
tion) that leads to the implicit Eq. (7), which is mani- 
festly finite and gauge and renormalization group invari- 
ant. The soft(er) scale in the coupling, too ~ y/aT <C T 
for a C 1, yields stronger screening than previously es- 
timated from (1). Remarkably, both the novel result 
(8) and the underlying running coupling (at T = 0) are 
consistent with corresponding 1QCD results. One might 
speculate whether Eq. (8) is especially suited for extrap- 
olation by its formal amenities - higher order terms could 
basically cancel each other if pcrturbative expansions 
have similar properties as asymptotic series [8]. 

Exemplifying the renormalization at T > is (at the 
considered level) particularly straightforward for the De- 
bye mass. Yet this quantity is a regulator for many ob- 
servables with an IR sensitivity and thus central in statis- 
tical field theory. Accordingly, there are numerous phe- 
nomenological implications for heavy ion physics. Men- 
tioned here explicitly with regard to the establishing no- 
tion of a sQGP is that the electric sector of QCD indeed 
seems to be quite strongly coupled, ole ~ a(mz>) ~ 0.6 
near T c , see Fig. 4. In the magnetic sector, on the other 
hand, the effective coupling apparently remains much 
smaller, «m ~ 0.3 near T c [9]. This indicates a rather 
nuanced picture of the sQGP, see e. g. [10]. 
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